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Abstract 



We consider reaction-diffusion systems on the infinite line that exhibit a family of spectrally stable spatially 
periodic wave trains uo(kx — Lot; k) that are parameterized by the wave number k. We prove stable diffusive 
mixing of the asymptotic states Mo {kx + 4>±;k) as x — > ±00 with different phases <j>- =fc 4>+ at infinity 
for solutions that initially converge to these states asi-> ±00. The proof is based on Bloch wave analysis, 
rcnormalization theory, and a rigorous decomposition of the perturbations of these wave solutions into a phase 
mode, which shows diffusive behavior, and an exponentially damped remainder. Depending on the dispersion 
relation, the asymptotic states mix linearly with a Gaussian profile at lowest order or with a nonsymmetric 
non-Gaussian profile given by Burgers equation, which is the amplitude equation of the diffusive modes in 
the case of a nontrivial dispersion relation. 
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1 Introduction 



We consider spatially extended pattern-forming systems that exhibit periodic travelling-wave solutions u(x, t) — 
uo(kx — ut; k) for a certain range of wave numbers k £ (ki, k r ). The profile Uo(9; k) is assumed to be 27r-periodic 
in 9 = kx — wt, where the wave number k and the temporal frequency uj are assumed to be related via a nonlinear 
dispersion relation to — uj{k). Examples are the Taylor vortices in the Taylor-Couette problem, roll solutions in 
convection problems, or periodic wave trains in reaction-diffusion systems. 

We are interested in the dynamics of perturbations of wave-train solutions of the above form. Since the lin- 
earization around a wave train always possesses essential spectrum up to the imaginary axis, we cannot expect 
exponential relaxation towards the original profile even for spectrally stable wave trains. Moreover, the periodic 
nature of the underlying wave train suggests that we should allow perturbations that change the phase or the 
wave number of the underlying profile. In these cases, we expect that diffusive decay or diffusive mixing of phases 
or wave numbers dominate the dynamics. In more detail, given a spatially periodic wave train ug(kQX — Woi; kg) 
we may consider (a) its diffusive stability, that is, its stability with respect to spatially localized perturbations, or 
else the diffusive mixing of the asymptotic states uo(k±x + </>±; k±) as x — > ±00 with (b) identical wave number 
fc_ = k + but different phases ^ cf> + or (c) different wave numbers fc_ ^ k + for initial data that converge to 
these states as x — > ±00. More precisely, consider an initial condition of the form 

u(x, 0) = u (q (x)x + c/)a(x);q (x)), q (x) ->• k ±1 cj> (x) ->• <f>± as x -> ±00, (1.1) 

where the functions qo(x) and <j)o(x) are bounded and small in an appropriate norm. We may then expect that 
the solution u(t, x) can, to leading order, be written in the form 

u(x, t) ss u (q(x 7 t)x + <j>(x, t) — uj t; q(x, t)), 

and the issue is to determine the behaviour of the phase <p(x,t) and the local wave number q(x,t) as t — > 00. 
As indicated above, we can distinguish three different classes of initial data, namely (a) constant wave number 
qo(x) = k and equal phases cj> + = 0_ at infinity for non-zero phase perturbations 4>o(x) ^ 0> which correspond 
to localized perturbations of the underlying wave train, (b) constant wave number qo(x) = kg but different phases 
cj) + 7^ 0_ at infinity, which correspond to a relative phase shift of the wave train at ±00, and (c) different wave 
numbers 7^ k + at infinity; see Figure [l] for an illustration. 

In this paper, we address the cases (a) and (b) for general reaction-diffusion systems 

d t u = Dd 2 x u + f(u) (1.2) 

with x £ R, t > 0, and u(x,t) £ K d , where D £ M. dxd is symmetric and positive definite, and / is smooth. We 
now outline our results and refer to Theorems [l] and [2] for the precise statements: 

(a) For localized perturbations of a single wave train, that is, for qo(x) = ko and </>_ = <f>+, we transfer 
existing stability results from specific systems Sch96 ( Sch98b, Sch98a. Ucc07j to general reaction-diffusion 
systems. In lowest order, the dynamics near a wave train can be described by the evolution of the local wave 
number q(t,x), and we prove that the renormahzed wave number difference ^[^(t 1 / 2 ^, t) — ko] converges 
towards a multiple of the x-derivative of the Gaussian ^J— - exp(— for an appropriate constant a > 0. 
This yields the asymptotics 

sup \u(x, t) — 0iim0* (x — c g t,t)dgu (6; k)\ < C2t~ 1+b as t — > 00, 

with <j)*(x,t) — yj=^ e_:E ^ 4at \ where 0ii, n £ K depends on the initial data, where a > and c g £ K 
are constants determined by the spectral properties of uq(-, ko), and where b > is a small, but arbitrary, 
correction coefficient. 
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(a) Diffusive stability 



(bl) Gaussian diffusive mixing of phases 




Figure 1: The panels illustrate different types of diffusive behaviour in the frame that moves with the speed of 
the group velocity c g . (a) Localized perturbations of wave trains decay diffusively like Gaussians; (bl) If u" = 0, 
then phase fronts develop when cfr- 7^ <f)+, and the wave number perturbation decays diffusively like a Gaussian; 
(b2) If u" 7^ 0, then phase fronts develop, and the wave number perturbation decays as determined by the Burgers 
equation; (c) Shown is the expected formal diffusive mixing of wave-number fronts in case uj" — 0. Solid and 
dashed lines indicate solutions at t = and for ( > 1, respectively, while the small solid curves in (bl)-(b2) 
indicate the amplitude- scaled spatially periodic wave train to visualize the phase shifts. 



(b) For perturbation that induce a global phase shift, that is, for qo(x) = ko and 0_ 7^ <fr+ but with 
\(/)d\ '■= \4>+ — <fi-\ small, we establish diffusive decay of wave-number perturbations. Specifically, the 
renormalized wave number converges to a Gaussian profile when w"(/co) — 0, while it converges to a 
nonsymmetric non-Gaussian profile when u/'(fco) 7^ 0. This latter case is the major result of this paper. 

The case (c) where qo(x) — > k± as x — > ±00 with fc_ 7^ k + is more difficult and depends crucially on the sign 
of uj"(ko). If a/' (fco) 7^ 0, diffusive mixing of the local wave number cannot be expected: instead, depending on 
the sign of u/'(fco)(fc + — k—), we expect that q(x, t) evolves either as a stable viscous shock or as an approximate 
rarefaction wave DSSS09 . If o/'(fco) = 0, nonlinear diffusive mixing can be expected, but, for some technical 
issues that we explain below, a rigorous proof remains open and is left for future research. 

The proof of diffusive mixing of phases of wave trains in systems with no S^-symmetry has resisted many 
attempts. With the rigorous separation of the phase variable (f> from remaining modes found in |DSSS09j . a 
new technique is now available to treat this question. This method combined with the renormalization group 
method |fjK92llFjKT94] . which has been applied for instance in |Sch96l IBch98bl lEgQOl iG^UOl I U ec041 IUec07l to 
a variety of pattern-forming and hydrodynamic systems, finally yields our results. Diffusive mixing results for 
the real Ginzburg-Landau equation, which has a natural decomposition into phase and amplitude variables due 
to its gauge symmetry, have been obtained for instance in BK92 , GM98J. 

The results in this paper were presented at the Snowbird meeting in 2007. Meanwhile, similar results on the 
diffusive stability of wave trains have been established in |JZ11[ IJNRZll] using pointwise estimates. 
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Notation. Throughout this paper, we denote many different constants that are independent of the Burgers 
parameters a,f3 and the rescaling parameter L > by the same symbol C. For m\,rri2 G N, we define the 
weighted spaces H m2 (mi)={u G L 2 (R) : ||«||fl- m a(mi)<oo} with norm \\u\\ Hm2{rni ) = \\up mi \\h™*(M), where 
p(x) = (1 + x 2 ) 1 / 2 and H m *(M.) is the Sobolev space of functions with weak derivatives up to order ni2 in 
L 2 (R). With an abuse of notation, we sometimes write ||tt(a;,t)||fl-m2( mi ) for the H m2 (mi)-norm of the function 
x i y u(x,t). The Fourier transform is denoted by JF so that u(k) := J r (u)(fc) = J e~ lkx u(x) dx for u G L 2 (R). 
Parseval's identity and J r (d x u)(k) — ifcu(fc) imply that T is an isomorphism between H m2 (mi) and H mi (m2), 
that is, the weight in physical space yields smoothness in Fourier space and vice versa. To indicate functions in 
Fourier space, we also write u G H mi {m2) instead of u G H mi {m2)- 



2 Statement of results 



2.1 Wave trains and their dispersion relations 



We assume that there are numbers ^ and uj G K such that (1.2 1 has a solution of the form u(x,t) = 
ua(kox — where uq{6) is 27r-periodic in its argument. Thus, uq is a 27r-periodic solution of the boundary- 
value problem 

k 2 Dd 2 g u + iod e u + f(u) = (2.1) 



with k = fc and cj = uq. Linearizing (2.1 1 at uq yields the linear operator 



Co = C(k ) = k 2 Q Ddl + LJ de + f(u Q (9)), 



(2.2) 



which is closed and densely defined on L 2 GI .(0, 2tt) with domain T>(Cq) — H 2 C1 .(0, 2ir). We assume that A = is 
a simple eigenvalue of Cq on L 2 er (0, 27r), so that its null space is one-dimensional and therefore spanned by the 
derivative dgUQ of the wave train. 



We may now vary the parameter k in (2.1) near k = ko and again seek 27r-periodic solutions of (2.1). The 



derivative of the boundary- value problem (2.1 ) with respect to u), evaluated at k = ko in the solution uq, is given 



by deuo- Since A = is a simple eigenvalue of Co on L 2 r (0,27r), we see that dguo does not lie in the range of 



£oi an d the linearization of the boundary-value problem (2.1| with respect to (u, uS) is therefore onto. Thus, 



exploiting the translation symmetry of (2.1) we can solve (2.1) uniquely, up to translations in 9, for (u,u>) as 



functions of k and obtain the wave trains 



i(x, t) — Uo{kx — uj(k)t; k), k G (fc;, k r ), 



(2.3) 



where w(fco) — w o and fc; < fco < fc r . In particular, wave trains exist for wave numbers k in an open interval 
centered around fc . We call the function fc H ► u)(k) the nonlinear dispersion relation and define the phase speed 
of the wave train with wave number fc by c p := w(fc)/fc and its group velocity by 



do; 
"dfc 



(fc). 



(2.4) 



To state our assumptions on the spectral stability of the wave train uq as a solution to the reaction-diffusion 



system (1.2), we consider the linearization 



d t v = C v 



(2.5) 



of (1.2) in the frame 9 = kox — u>ot that moves with the phase speed c p = ujo/ko- Particular solutions to this 

-wave ansatz 

u(0,t) =e x ^ t+iee / k °v{9,e), (2.6) 



problem can be found through the Bloch-wave ansatz 
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where £ £ R and v(0, £) is 27r-periodic in for each I. In fact, since v(d, £ + fco) = e 11? u($, £), we can restrict £ to 



the interval [— fco/2, ko/2). Substituting (2.6) into (2.5), we obtain 

C(£)v = X(£)v (2.7) 

with a family of operators £ given by 

C(£)v = e-" fl / fc °£(e" 9 / fc °5(0, £)) = fcgD (9 e + i£/k ) 2 i + cu(dg + i£/k ) v + f'(u (6))v, (2.8) 

each of which is a closed operator on Lp er (0, 2ir) with dense domain Hp C1 .(0, 2ir). In particular, £(£) has compact 
resolvent, and its spectrum is therefore discrete. We can label the eigenvalues of C{£) by indices j £ N and write 
them as continuous functions Xj(£) of £. In addition, we can order these eigenvalues so that Re A J+ i(0) < Re Aj(0) 
for all j. In fact, the curves £ H> Xj(£) are analytic except possibly near a discrete set of values of £ where the 
values of two or more curves Xj(£) for different indices j coincide. 

Next, we assume that Ai(0) is the rightmost element in the spectrum for £ = 0. Since we assumed that A = is 
algebraically simple as an eigenvalue of C, there is a curve Xi(£) of eigenvalues with Ai(0) = 0, and this curve is 
analytic in £ for £ close to zero. We call the curve Xi(£) the linear dispersion relation and denote the associated 
eigenfunctions of £(£) by v\(9,£). We shall compute the derivative dAi/ d£ and recover the group velocity as 
defined via the nonlinear dispersion relation, namely 

-Im^Ailfco = -c p + dkuj(k ) = -c p + c g . (2.9) 

We remark that the phase velocity c p appears in this formula solely because we computed Ai in the frame moving 
with speed c p , while u> was computed in the steady frame. We also note that the signs of the second derivatives 
of Ai and uj are, in general, not related. Finally, we assume that ReA"(0) < and that all other eigenvalues 
Xj(£) satisfy HeXj(£) < — ctq- The following hypothesis summarizes the assumptions we made so far. 



Hypothesis 2.1 (Existence of spectrally stable wave trains) Equation (1.2) admits a spectrally stable wave 



train solution u{x,t) = uq(8) with 9 = k$x — Lo$t for appropriate numbers fcg / and ojq £ K, where uq is 2tt- 



periodic. Spectral stability entails the following properties. First, the linearization £ of {1.2) about u has a 
simple eigenvalue at X — 0. Furthermore, the linear dispersion relation Xi(£) with Ai(0) = is dissipative so that 
A"(0) < 0, and there exist constants ao,£o,ao > such that KeX\(£) < — cro for \£\ > £q and KeX\(£) < —olq£ 2 
for \£\ < £o, while all other eigenvalues Xj(£) with j > 2 have Re Xj{£) < —<jq for all £ £ [— fc/2, fc/2). 

Standard perturbation theory yields that the wave trains u (kx — u{k)t; k) are also spectrally stable, possibly 
for a smaller interval fc; < fc; < k < k r < k r of wave numbers than the interval of existence. By changing ki, k r 
accordingly, we shall assume from now on that the wave trains uo(-;k) with k £ (ki,k r ) are spectrally stable 
with uniform constants £o,o~o,oto. 

For later use, we collect a few properties of the linear dispersion relation and refer to [DSSS09, §4.2] for their 
derivation. We denote by 

£ ad w = k^DdgU - u d e u + f'(u (8)) T u 
the Lp Cr ((0, 27r))-adjoint of Cq = £(fco) and let w a <i be a nontrivial function in its null space with the normalization 

(M a d,5 e W >L2(o,2 7 r) = 1- (2-10) 

Using the adjoint eigenfunction, we have 

A;(0) = i(ii ad ,c p a 9 wo + 2k Dd^u ) L2 - i(c p - c g ) £ iM, (2.11) 
Ai'(0) = - (u ad , Ak Q Dd k d e u + 2Ddgu ) L2 £ R. (2.12) 

We shall also use the identity 

dtv 1 {;0)=id k u (2.13) 

that was established in [DSSS09, §4.2]. 
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2.2 Statement of results 



Throughout this section, we fix the wave number ko of a wave train uo(kox — Lo^t;ka) of the reaction-diffusion 



system (1.2 1 that satisfies Hypothesis 2.1 We then set 



w =w(fc ), Cp^uj /k , c g =uj / (k ) 7 P = --u)"(k ), 9 = k x-uJot 

and write 

\ 1 (i) = i(c p -Cg)£-af +0(e 3 ) (2.14) 

for the expansion of the linear dispersion relation of uo(-;ko). For convenience henceforth we write k = fco- 
Before we state our result, we remark that the decomposition of the initial data that we shall use below in the 
statements of our theorems is not unique. This non-uniqueness will be removed in the proofs but does not affect 
the conclusions made in the results below. 

Our first result states that uq is diffusively stable with respect to localized perturbations and extracts the leading- 
order behaviour of the displacement for large times. For notational convenience, in the following we consider 
initial conditions at t = 1. 



Theorem 1 (Diffusive stability) Let Uq(-; k) be a spectrally stable wave train that satisfies Hypothesis 2.1 and 
pick b € (0, 1/2); then there are e,C > such that the following holds. If, for some 8q <e [0,27r), 

u{x,t)\ t= i = u a (6 - O + (j> (x); k) + v (x) with \\<f>o\\m(3)i IN||iP( 3 ) < e, (2.15) 

then the solution u(x,t) of \1.2ty exists for all times t > 1, it can be written as 

u(x, t) = u (9 — Q + <f>(x, t);k) + v(x, t), 

and there is a constant <j)\i m € K depending only on the initial condition so that 

supU(a;,t) - <hmG(x - c g (t - l),t)\ + \v(x,t)\ < Cr 1+b , (2.16) 

where 

G(x,t) = —^e- x2/(4at \ (2.17) 



In particular, we have 



sup \u(x,t) - u (6 - 9 + 4>n m G{x - c g (t - l),t);k)\ < dt 



-i+b 



Next, we discuss diffusive mixing of phases for non-localized phase perturbations. In this situation, the precise 
asymptotics of perturbations depends on [3 = —^oj"(k). 

Theorem 2 (Diffusive mixing of phases) Let uq(-; k) be a spectrally stable wave train that satisfies Hypoth- 
esis 



2.1 and pick b € (0, 1/2); then there are constants e, C > such that the following holds. 



(i) Assume that (3 = —^(jj"(k) = and u(x, t)|t=i = u o (0 — 6 + (j>o(x); k) + Vq(x) with 4>o(x) — > <j>± for x ±oo, 
\4>d\ = \4>+—<j>-\< s, and 

Uo(-)\\HH2),\\vo\\H H 2)<e. (2.18) 
Then the solution u(x,t) to \1.2\ exists for all t > 1, and can be written as 

u(x,t) = u (8 — &o + <j>*{x — c g (t — 1),<); k) + v(x,t), 
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where 



<f>*(x,t) = ^ 



(p-)evi (x/v at), with erf(x) 



1 

In 



(2.19) 



and sup xm \v{x,t)\ < C t'^+K 



(ii) The same result holds if ft = —^uj"{k) ^ 0, with <j)*(x,t) replaced by 



(jf[x,t) = - ln(l + zeri(x/Vat)), ln(l + z) 



(2.20) 



Remark 2.2 Clearly, the decompositions (2.15) and (2.18) are not unique. For instance, 0o = would be 



one possibility in (2.15), but we may shift perturbations between 0o an d fo- In the proof we shall fix this 



non-uniqueness via mode-filters. 

The higher weight in the initial conditions in Theorem [l] vs. [2] is due to the fact that in Theorem [I] we want 
to extract higher order asymptotics, i.e., faster decay. The asymptotic phase-profiles in (2.17l,( 2.19[ ) and (2.20) 



only depend on k via a from (2.14) (and on /3 for (2.20)). In particular, they are independent of the phase-speed 
c p and therefore are formulated in x and t. 

Remark 2.3 Formally, we may as well describe the diffusive mixing of wave numbers in case to" — 0, see 
Remark 2.7 However, then the rigorous separation of the (then unbounded) phase, see ^j3j becomes more 
difficult. Therefore, we will not consider this case here. 

2.3 The idea 

The translation invariance of (|1.2[) and the fact that by assumption we have periodic wave trains uo(6;k) for 



wave numbers k in a whole interval (fe;,fc r ) suggest to consider initial conditions for (1.2) of the form 

\uic(x) — uo(k±x + 4>±] k±)\ — > as x — > ±oo. 



(2.21) 



The behavior of the corresponding solutions can be discussed formally if we assume that the initial phase shift 
4> + — <fi- or the initial wave number shift q + — g_ happens on a long spatial scale. We make the ansatz 

u(x, t) = u (k Q x - LJ t + $(X , T);k + 5d x <S>{X, T)) (2.22) 

where < S <Si 1 is a small perturbation parameter that determines the length scale over which the wave number 



is modulated by the function dx&, and where X and T are long spatial and temporal scales. Plugging (2.22) 



into (1.2) and comparing equal powers in 8 it turns out (see |DSSS09l §4.3]) that 

(X, T) = (5(x - c g t), 8 2 t) , where c g = w'(fc), (2.23) 
are the correct spatial and temporal scales, and that q(X,T) := dx&(X, T) should satisfy the Burgers equation 

d T q = ad 2 x q + Pd x (q 2 ), a = -^£(0), = -\<J'{k), (2.24) 
while the phase &(X, T) itself satisfies the integrated Burgers equation 

= ad x <f> + /3(dx$>) 2 - (2.25) 
Note again that in case 13 = (Theorem [2] (i)) ( pT24| resp. (I2.25I) are the linear diffusion equation. 



Two questions arise: (a) What do we (formally) learn from (2.24) resp. (2.25)? (b) In what sense, i.e. in what 



spaces and over what timescales, do solutions of (2.24) via (2.22) approximate solutions of (1.2), and can we give 
rigorous proofs for that? 
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To answer (a) we briefly review some well-known results about dynamics and stability in the Burgers equation in 
the following section. With this in mind we turn to (b). One way to translate the formal analysis into rigorous 
results is to give estimates for the difference between the formal approximation 

U apptm (x, t) = u (9 + *(X, T;8);k + 5d x $(X, T; 6)) 



and a true solution u(x,t) of (1.2), on sufficiently long time scales. In [DSSS09 this has been achieved for 
a variety of cases using a separation of the critical mode (the phase mode) from the exponentially damped 
remaining modes by Bloch wave analysis. Here, for special initial data we obtain the diffusive stability results 
and the mixing results from Theorems [I] and [2j The proofs heavily rely on the coordinates from |DSSS09j . which 
are introduced in Sj3l 



2.4 Dynamics in the perturbed Burgers equation 



In the spectrally stable case the amplitude equation for long wave modulations of the local wave number is given 
by the Burgers equation. For given classes of initial conditions the behavior of solutions of the Burgers equation 
is well understood, and, moreover, this behavior is stable under perturbations of the Burgers equation. Thus, 



before proving our results for (1.2) with initial data (2.21) we briefly review some well-known results about 



the (perturbed) Burgers equation, cf. [BKL94. MSUOlj and |Uec07[ Sec. 3]. This also motivates the ideas and 



methods of the proof. To keep track of a and /3 we do not rescale (2.24) to the standard form d T q = d?q + d^(q 2 ). 



c solves drv = ad x v 



The Burgers equation (2.24) has Galilean invariance: if q solves (2.24), then v — 
2cf3dxv + f3dx{v) 2 which can be transformed back to (2.24) via Ih>I + 2cf3T. Thus, concerning the stability 



of constant solutions of (2.24) we can restrict to q = 



We add a higher order perturbation in the form of a total derivative to (2.24) and for notational convenience we 
take initial conditions at time T = 1. Thus we consider 



d T q = ad x q + f3d x {q 2 ) + jdxh(q,d x q), q\r=i = Qo, 



(2.26) 



where for simplicity h(a,b) = a dl b d2 is a monomial. For 7 = we again have the Burgers equation. The 
perturbation is assumed to be of higher order. To make this precise we define the degree 



dh = d\ + 2c?2 — 3, and assume that di < 1 and dh > 0. 



(2.27) 



The mean J R q(X, T) AX is conserved also by the perturbed Burgers equation (2.26 1 . Diffusive stability of q = 



in ( 2.26[ ) is based on the fact that solutions to the linear diffusion equation in Fourier space concentrate at wave 
number k = 0. Roughly speaking, for initial data in i 1 (IR) that decay like |^| _rl , the solutions of dxq 
fulfill 



ad 2 x q 



q{X,T) = T-V+V^q^^HjiX/VT) + 0{T- n ' 2 ) for T -> 00, 



(2.28) 



3=0 



where Hj is a multiple of the (scaled) j Hcrmitc function H(x) = (— l) J d%, cxp(— x 2 /(4q)). Thus, if <7o(0) = 
± J R qo{X) dX £ then ||g(.,T)|| L « < CT-^\\q Q \\ L ^ while for q(Q) = we have \\q(; T)|| L » < CT~ 



\QoWl 1 



In the second case it turns out that solutions to the nonlinear equation (2.26) with zero mean have the 
same asymptotics as solutions to the linearization with zero mean. Thus, both nonlinear terms /3dx{q 2 ) and 
jdxh(q, dxq) are called asymptotically irrelevant. 



For q(0) 7^ only jdxh(q, dxq) is irrelevant, and there is a nonlinear correction to the dynamics for (2.26) 



compared to (2.28). To derive this we use the Cole-Hopf transformation 

^ y/a d Y Q(Y, T) 

P Q{Y,T) ' 



Q{X,T) =exp 



q(Y,T) dY 



q(X,T) 



Y = Xjsfa, 
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which transforms (2.26) with 7 = into the linear heat equation drQ = d x Q, Q\t=i — Qo, with limx-»-oo Qo(X) 
1 and limx^oo QoPQ = 1 + 2 > 0, i.e. 



ln(l + z) 



q(Y, 1) dF. 



Since limT->-oo Q (^\/TX,Tj = l + zeri(X) + 0(l/VT) we find that the solution q to the Burgers equation (2.261 
with 7 = satisfies 



T lim Vfq (VfX, T) = ln(l + z evf{X / yfa)) =: f* z (X), (2.29) 

with rate 0(1/Vt). Therefore, if (3 7^ 0, then the renormalized solutions converge toward a non-Gaussian limit 



f*(X). Again, the same behavior can be shown for (2.26) with 7 7^ 0. We summarize these results as follows 



Proposition 2.4 For each b € (0,1/2), there exist C\,C 2 ,Tq > such that for solutions q of the perturbed 



Burgers equation (2.26) the following holds 



i) Assume that \\qo\\H 2 (i) < Ci and Qo(X) dX = 0. Then there exists a qn m G R such that 



Tq (VfX,T) -q linl Xe- x2 / 4a 



< C2T-2- 



ff 2 (3) 

Thus, \\q(X,T)\\ L i < C 2 T-^ 2+b and \\q(X,T)\\ L ~ < C 2 T- 1+b . 

ii) Assume that A = J_ qo(X) dX 7^ 0, j3 — 0, and ||<Zo||h 2 (2) < Ci- The 



(2.30) 



T 1/2 q (VTX,t) - 



A 



\J Ana 



-X 2 /(4a) 



< C 2 T-^ +b , 



iP(2) 



(2.31) 



and consequently \\q(X, T)\\l°° < C 2 T x l 2 . 

Hi) Assume that A = qo(X) dX 7^ 0, j3 7^ 0, and ||9o||_ff 2 (2) < Ci- Then 

T i /2q ( Tl/2x T \ f {x) C2T -x/2 +6) where r {x) = 

\ J ff 2 (2) [3\jAn 1 + zcrf {X/y/a) 

id ln(l + z) = § q (Y) dY. In particular, again \\q(X, T)\\ L =c < C 2 T" 1/2 . 



(2.32) 



Remark 2.5 a) By translation invariance of (2.26), we can replace go in Proposition 2.4 by qo(- — Xq) for some 



Xq £ R and obtain the corresponding results for q{X — Xq,T); w.l.o.g. we set Xq = 0. 

b) The higher weight for qo in Proposition |2.4[ i) compared to (ii),(iii) is due to the fact that we want to isolate 
higher order asymptotics (with faster decay). For this we need F(qo) £ H 3 (2) C 2 . 

c) The profiles in (ii), (iii) are explicitly given in terms of A due to the conservation of J qdx, i.e., since the 



right-hand side of (2.26) is a total derivative. On the other hand, the constant qn m in (i) in general depends on 



qo in a complicated way. 

d) The local phase $, which is related to the wave number q by q = d x $, satisfies the (perturbed) integrated 
Burgers equation 

dT<f> = ad 2 x <S> + (3(dx<S>) 2 +lh(d x t>,d 2 x <f>), 1) = $ {X). (2.33) 



For (2.33) there exist C\, C 2 > such that we have the following asymptotics. 
i) If 1 1 $0 1 1 h 3 (3) < Ci then there exists an (f>n m = — 2aqu m € R such that 



T 1/2 $ (VrX,T 



5ii m e 



-X 2 /4a 



< C 2 T- x l 2+h . 



(2.34) 
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Thus, the renormalized phase converges toward a Gaussian. 

ii) If /? = and $ a (X) -> $± as X -> ±00 with |$ + - $_| < C x and ||$o|| ff =(2) < C u then 

||$(Vf X,T) - $*(X)|| ff 3 (2) < C 2 T- 1/2+b 
where = $_ + ($+ - $_)erf(X/ 1 /a). 

iii) If /3 ^ and and $o(-^) -) $± as I ^ ±00 with |$ + - $_| < Ci and || < I ) ol|if 2 (2) < Ci, then 

\mVfx,T) - $*(x)|| ff3(2) < c 2 r~ 1 ^ 2+b 

where = $_ + | ln(l + zexi(X/ ^/a)), ln(l + z) = £($+ - $_). 



J 



Remark 2.6 We briefly want to explain the reason for (2.27) and the idea of (discrete) renormalization. If 

(2.35) 



/ qo(X) dX 0, then, for L > 1 chosen sufficiently large, we let 

q n (Z,T)=L n q(L n li,L 2n T). 
Then q n satisfies 

d T q n = adf <?„ + Pd^ql) + 'fL^ n d i h n (q n , d ( q n ), 



with 



where dh 



(2.36) 



(2.37) 



h n (q n ,d ( q n ) = L~^ +2d >-V n q d n i(d ( q n ) d \ 
2d 2 - 3 > due to (|2.27|). Next, solving d T q = ad 2 x q + (3d x {v 2 ) + jd x h(q, d x q) for T £ [1, 00) 



is equivalent to iterating the renormalization process 

solve (|2.36| for r € [L~ 2 , 1] with initial data g„(£, L~ 2 ) = L<2„_i(L£, 1) € A", 



(2.38) 



where A" is a suitable Banach space. Since (formally) L^ n d^h n in (2.36) goes to zero, in the limit n — > 00 we 
recover the linear diffusion equation (if f3 = 0) respectively the Burgers equation (if ft ^ 0) for q„, with the 
known asymptotics (2.28) respectively (2.29). Similarly, if J qodX = 0, then we scale 

q n (Z,T)=L 2n q{L n Z,L 2n T), (2.39) 

and (independent of whether f3 is zero or not) end up with the linear diffusion equation in the respective 
renormalization process. To make this rigorous we need a suitable Banach spaces X and rigorous control of the 
iterative process (2.38), and again we refer to |BKL94| and |Uec07[ sec. 3] for details. However, two observations 
are most important, (a) In (2.37) we see that each derivative in x gives an additional L^ 1 in the rescaling. (b) 
The diffusive spreading in physical space corresponds to concentration at n — in Fourier space according to 
T(Lu(L-))(k) = u(k/L). Thus, for the linear part, only the parabolic shape of the spectrum A(k) = —an 2 of 
ad 2 near k = is relevant. 

Remark 2.7 If q (X) ->• q ± for X ->■ ±00, then q(VTX,T) = Q* a {X) + C(T -1 / 2 ) as T ->■ 00 for the solutions 
of <7t = adxq, where Qq(X) = g_ + (q + — q_)eit{X / ^fa) . Thus we have diffusive mixing of the wave numbers. 
Then, for /? = and for suitable qo, we have the asymptotics 



q(VfX,T) = Q*{X) +0(T~ 



1/2N 



as T 



DC 



(2.40) 



for (2.26), where \Q* (X) — Qq(X)\ < CeT x / 4 , i.e., we have essentially the same asymptotics as in the linear case, 
with a small localized correction, see [BKL94] . On the other hand, for (3 7^ a front is created, see UDSSS09 . 
However, here we do not further comment on this case since below we focus on diffusive mixing of phases. 
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3 The separation of the wave numbers 



3.1 The ansatz 



Only special systems such as the cGL have an 1S 1 -symmetry and therefore a natural decomposition into amplitude 
and phase. Hence, the first step is to extract from a general reaction-diffusion system an equation for the phase, 
and then out of this for the wave number. We follow the formal derivation made in DSSS09 which uses a 
multi-scale expansion which however we cannot assume a priori. Thus, here we proceed as follows for the 
reaction-diffusion system (1.2). As above we change coordinates via — kx — cot, and obtain 

d t u = k 2 DdgU + udgu + /(it). 



(3-1) 



(3.2) 



(3.3) 



A stationary wave train uo(9; k) of (3.1) with period 2tt satisfies 

k 2 Ddgu Q + ud e u Q + f(u ) = 0. 
Given a smooth phase function <j)(d,t) we seek solutions of the form 

u{9, t) = u (&, k(l + fl^(0, *))) + u>0?, *). 
where the phase (/>(■&, t) and the coordinates and ■& are related by 

= i? t). (3.4) 

Roughly speaking we require that d®4> is small, uniformly in and that 0($, t) is close to the asymptotic profile 
we want to extract. Still, (3.3) adds an additional degree of freedom by introducing <f>] we later add additional 
conditions on <f> and w, via mode filters, to remove this additional degree of freedom again. 

Remark 3.1 It might seem more natural to make the ansatz 

u (6, t) =u (9 + <j){6, t);k(l + d e <p(6, t))) + w{6, t) 



(3.5) 



instead of (3.3 1 . However, we need to be able to relate the dynamics of u(9, t) back to properties of the wave-train 
«o(#; k). Thus, we would need to express u(0, t) in terms of & = 9 + (f>(9, t), i.e., 

u (9 + 0(0, f); fe(l + dg<j)(9, t))) i— >■ u ($; fc(l + dgcj>{6(§, i),t))) 



which involves the inverse 9{d, t) of the function d = 9 - 
the forthcoming analysis much more complicated. 



, t). The occurrence of this inverse would have made 

J 



Remark 3.2 Suppose that we found a phase function 4>{d, t) with small derivative 9#0($, t) so that ( 3.3 ) satisfies 
(3.1). Using the implicit function theorem, we can then, a posteriori, solve (3.4) for d as a function of 9 which 
(9,t), where 



M),t). 



,t),ty,k(l + dg<j>(6 + 4>(6,t),t))) 



is of the form ■& = 9 + 

4>(0,t) = ct>(#,t) = < 

In particular, we see that 

uo(0;*(l + 080(iM)))=«o(0 + ^(0 + ' 
and 

< + 0(9,t),t) = (l + dg]>(d,t))d <t>(d + 4>(9,t),t) 

- d g ^(9 + fo9, t),t) + O(\d e 0(9 + 0(0, t), t)\ 2 ). 

, t) replaced by 



(3.6) 



d 

d0' 



Thus, to leading order, the solution (3.6) is of the desired form (3.5) with 



ht),t). J 
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We now substitute the ansatz (|3.3[) into (3.1) and derive the resulting PDE in <&. We use the notation 



ug:=tio(0;fc(l + a^)). d j 4:=(d j u )'t' :=(d jUo )(&,k(l + d^)), j = 4,k. (3.7) 

Assuming that d§4> is small, we obtain 

dd 1 di? -dt<t> 



and therefore 



and 



Thus 



d6 l-dfitf dt l-d#<j>' 

d Id d 2 / 1 d x 2 



dO l-d^dti' d9 2 \l-d#4> dti 



d 2 u fid kdl4> d x 1 



d6 2 \l-d#<t> dtf l-d#<j> dkj °' 

du _ -dt<f) ^ / dl4>d t A 



dt l-d#6 \ 1-d* 



d#ul + k — — — + dodt4> d k ul 



dw dw dw dt(t> dw 1 dw d 2 w ( 1 d N 



dt dt d-dl-d^ dG l-d^dti' dO 2 \l-d#<j)d<d 



w. 



d$u% - k I 8%<f> - 9fl9 t ) d k ul + d t w - d^w (3.8) 



i 9 fc<9 2 <?(> a\ 2 , / i a \ 2 



2 n x " , ^-dS 



k z D 



l-drfdd 1-d^dkJ u \l-dt4>d0 
-{k?Dd%u Q - u>d#u + f(u Q )) + f(4 + w) 



where we used (3.2 1 in the last equation. 



Our goal is to separate the critical modes, which involve the dynamics of <j), from the damped noncritical modes 
using the eigenfunctions of the linearization C(k). This is done via Bloch waves which we introduce next. 



3.2 Bloch wave analysis 

Bloch wave transform J is a generalization of Fourier transform J- . We briefly review the main properties and 
refer to |RS78| ISch98b[ ISca99 , DSSS09 for proofs and further details. From now on we use a slightly rescaled 
Fourier transform, namely 

-i /*oo /*oo 

wit) = (Fw)U) = — / e-' m ' k wm dtf, wCff) = (T^w)^) = / e M/k w(£)d£, (3.9) 
2nk J_ 00 

and thus, denoting the classical Fourier transform (where k = 1) by .7-1, 

[Fw){l) = \{F x w){l/k) and (Ju>)(i?) = {T^w){§/k). (3.10) 
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Then, for sufficiently smooth and rapidly enough decaying functions w, we have 

-fc/2 

^ t+jU) l k w{t + jk) d£ 

r k/2 



(J- 1 ™)^) := w($) = f e iM/k w(l) dl = ^ / 

J — OO „■ f- 17 J — 



iM/k 



k/2 



dl = 



k/2 



-k/2 

m/k w($j) dl 



(3.11) 



where 



w(Q, I) = (Jw)(&, t)=Y^ e ij{> w(e + jk). (3.12) 

iez 

Similar to the Fourier transform, the Bloch transform can be defined for tempered distributions. By construction, 

u)(tf + 27r,£) =w{ti,l) and w{d, I + k) = e^w(§, I), (3.13) 

such that we can restrict ourselves to I £ [—k/2, k/2). The Bloch transform of the product of two functions Wi 
and w 2 in $-space is given by the convolution 



J[wfW 2 }(#,£) = kDi *u> 2 ](iM) = 



k/2 



k/2 



(3.14) 



of their Bloch transforms W\ and w 2 in Bloch space, where (3.131 is used for \t—t\ > k/2. The analytic properties 

ty 

u(d,l)\ 2 dl dl?. 



of the Bloch transform are based on a generalization of Parseval's identity 

r 2ix c-k/2 



\u{$)\ 1 dl? = 2irk 

lQ J-k/2 

As a consequence, Bloch wave transform is an isomorphism between H m2 (mi), and the space B mi (m 2 ) of 
functions u(d,l) that are 27r-periodic w.r.t. i?, satisfy (3.13), and whose norm 



mi rn 2 r 2ir fk/2 



\ u \\B m i(m 2 ) 



EE 



\didiu{ti,l)\ 2 dl di? 



j=0i =o Jo J - k / 2 

is finite. We now collect a few more properties; see, e.g., |Sch98bj or |DSSS09l §5.2] for more details and proofs. 

Remark 3.3 a) If iui(i?) is 27r-periodic in i? and the support of the Fourier transform w 2 of a complex- valued 
function 102(1?) lies in (—1/2, 1/2), then we have 



J[w 1 w 2 ]{$,l)=w 1 {ti)w 2 {l). 



(3.15) 



Due to (3.15 1 Bloch transform is useful to analyze differential operators with spatially-periodic coefficients, which 
in Bloch space become multiplication operators. 

b) Since we are interested in functions which do not necessarily decay to zero at infinity, we employ a method 
already used in [Sell!) 1 to extend multiplication operators from the space L 2 of square-integrable functions to the 

space L 2 ^ of uniformly locally square-integrable functions equipped with the norm ||w|j L 2 = sup / |m(?/)| 2 dy. 

ul x£R Jx 

We recall that 

H% = {u : R -> M; \\u\\h™ = \\u\\ H ™(x,x+i) < °o with lim ||u - T v u\\ H % -> 0} 
where [T y u](x) = u(x + y). Now let m, s € Z with m + s > and m > 0, and consider a function 

M: M^i(i/"+ s (0,2^),i7™ r (0,27r)), I ^ M{1) 
which is C 2 in the Bloch wave number I. Then M defines a bounded operator M. : H"\ +s — > with 

\\ M \\l(HZ + \H^) < C ( m ^)W\\cl((-k/2.k/2)MH^ s ,H ?CI ))- ( 3 - 16 ) 

Clearly, this can be extended to multi-linear operators. 
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3.3 Mode filters, and separation into critical and noncritical modes 



Our goal is to separate the dynamics of the eigenmodes ?)i(i9, £) associated with the critical eigenvalues \i(£) 
from the remaining modes, which are linearly exponentially damped and therefore called non-critical. We use 
mode filters to obtain this splitting. 



Due to Hypotheses 2.1 there exists a number £\ with < £\ <C 1 so that the eigenvalue \\{€) of C{£) is bounded 
away from the rest of the spectrum for \£\ < l\. Therefore, there exists an £(£)-invariant projection 



Q c (£) 



27ri 



[A - C^)]- 1 dX 



onto the space spanned by Vi(i!),£), where T C C is a small circle that surrounds \\(£) counter-clockwise in the 
complex plane and does not intersect the rest of the spectrum of £(£) for this fixed £. For £ = we have 

Q c (0)£(-,0) = («ad,«)«i(-,0), 

and similarly Q c (£) can be expressed by using the scalar product with the adjoint U a d{',£) in Bloch space. 
We choose a nonincreasing (for £ > 0) C^°-cutoff function % : R — >• [0, 1] with 



x(£) = 



1 for \£\ < 1, 
for \£\ > 2, 



(3.17) 



and define 



Pg(£) = Q c (£) x 
P^(£) = Q°(£) x 



and 



p c {£) = Q c (£)x 

These operators commute for each fixed £ and satisfy 



F%(£) :=1-Q C (£) X 
P m fW :=l-Q c (£) X 

P s (£) :=l-Q c (£) X 



(1 - P<)Pg = = (1 - Pg)P^ t , (1 - P s )P f l = = (1 - P S )P^, 

P c -4- P s — 1 P c 4- P s — 1 

We define scalar- valued operators p%{£) and p^{£) implicitly by 



(3.18) 



(3.19) 



3.3 d) implies that each of the operators above extends to a bounded operator on H"\ S . The resulting 



Remark 

operators will be denoted by the same letter but with the superscript ~ being dropped. 



The mode filters and Pl^ are now used to separate the critical and noncritical modes in (3.8), while Pf s and 
Pf s are used to limit the Fourier support of the critical modes. We write (3.8) given by 



d t <P 
1 - d. a q> 

= k 2 D 



d-aul -k[d, 



1 - d#(f> 



d$d t 4> d k % + d t w 



d kd#<f> d 



1 - d#(t> d$ 1 - d$<j) dk 

w - — — — (<9^Mq + k(dl<l>)d k u$ + da 
l — o$<p \ 

(k 2 Dd 2 u Q - utdtfUo + f(uo)) + f(v$ + w 



1 



9t4> 
1 - d#( 



d 



1 - d*6 dd 



(3.20) 
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MS 



where 



B d t <t> 
B 1 (dtf(/),w)dt(i> 



[—B Q + Bi(d^4>,w)]d t (t> + d t w = -C t d§4> + £w + G(d$<t>,w), 
(d$u Q - kd k u dv)d t (t>, 

—£(kd$(j)dkUo) + k 2 D(2d$(f>dluo + d^d^uo) — cod^(j)d^uo 
k [£(id$4> d t vi) + kD{2d$4>dlu + d^cj)d^uo) - c v d$<l> d$uo\ , 



(3.21) 



(3.22) 



d$u 



1 - doc, 



1 - / 1 - 



and where G is contains the remaining terms. In the calculation above, we used that dgVi = ic\.ito, see (2.13). 
The symbol Ci is used since in the critical modes d§Ci corresponds to C, see (3.33) below, i.e., Ci resembles an 
integration of C. Clearly, 

Bt(d ct>,w) = O(\d^\ + \w\), G(d^w) = (D(\d^\ 2 + \w\ 2 ). 



Our goal is to replace (3.21 ) with the system 



d t PgB <t> = PgCM + P^ { Bi{da<t>, w)d t <t> - P^G(d^ w), 

d t w = Cw + P! s B d t (j ) -PlL l d^^P s m{ B l {d^w)d t ^ + P i mi G{d^,w) 



(3.23) 



for [4>,w). Subtracting the first from the second equation and using ( |3.18[ ), we see that solutions of (3.23) give 

(3.24) 



solutions of (3.21). Alternatively, we may consider the system 

p%dd$(j) + p c mi Bi (d#<p, w)d t 4> - p c mi G(d$4>, w), 
Cw + P f s s B d t (b - PlUd^ - P^Btid^, w)d t 4> + P s mf G{dv4>, w 



dtPfsBofj) 
d t w 



for (cj),w), where the first equation is now scalar- valued. Inspecting (3.19) we see that (3.23) and (3.24) are 
equivalent. We shall require that (</>, w) satisfy 



and 



(1 - P s )w = 



(3.25) 



(3.26) 



for all t > 1. Since P s commutes with C, it follows from (3.18) and (3.24) that (3.26) holds for all i > 1 if it is 
true for t = 1. 



It remains to check whether (3.25) is respected by (3.24) and to calculate the operator p c ls Bo to see whether 
(3.24) is a proper evolution equation. Due to the properties of the multiplier p c mi , we know that 

suppF [p c ni( (B 1 (d$<p,w)dt<f> - G(d$ct),w))} <e X 



for any sufficiently smooth function (/). ^From (3.22) we find that the operators Bo and % have 27r-periodic 
coefficients in $ and are multipliers in Bloch space which allows us to use Remark 3.3 For any function <fi that 
satisfies (3.25), we then obtain 

PgJ [B <t>] = Ff a (£)J[B <t>}(#,t) O ${£) X (U/h) Q c (£) (d#u (#) + 0{i)) 



= mx{^/ti)^ + 0{t))v 1 (d,£)^ (1 + 0(h)) <!>{£) v x {$,l) 



where the 0(fi)-term is a multiplier and (1 + 0(l{f) < 



has support in X. Therefore, using the definition (3.19 1 



of Pf s and denoting the operator associated with the 0(£i)-term by B 2 , we get 

P C M= (1 + B 2 )4> 



(3.27) 
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for all 4> that satisfy (3.25), where B 2 has norm H-B2II = 0(£i) and respects (3.251, i.e. suppj r [i?20] C I. Since 



similar arguments apply to the multiplier Ci, (3.25) is indeed preserved by (3.24) 



For all (4>,w) for which (d$4>,w) is small and <fi satisfies (3.25), the first equation of (3.24) can be written as 



d t <j> = [1 + B 2 + p^Bxidtfcf), w)} 1 b4A<9,j0 - p^G(d^, w)} . 

Substituting this expression for dt4> into the second equation of ( |3.24[ ) for w, we arrive at the system 

d t 4> = [1 + B 2 + p^B^d^, w)]- 1 \p%Cidi><j> - P c m{ G{d^ to)] , (3.28) 
d t w =Cw - Pgdd^ + P^G(d^, w) 

+ [PgBo-P^Btid^, w)} [l+B 2 +p c mi B 1 (d^, w)]- 1 [p c ts C t d^-p c mf G{d^ w)} . (3.29) 

Thus we have a splitting of the critical modes (f> and the noncritical modes w. 

3.4 The system for wave numbers and damped modes 

We now replace </> by ip = d$4> and obtain 

W =di> [1 + B 2 + P c m{ B 1 (^,w)}- 1 \j^Ci^-j^tG{il>,w)] , (3.30) 
d t w =£w - Pldii + P n s lf G(^, w) 

+ [PgBo-P^Bxi^w)] [l+B^p^B^.w)]- 1 \pldiJ-p c m{ G(i>,w)] , (3.31) 

which we also write in short as 

d t V = AV + F{V), (3.32) 
where V = (ip, w), A is a linear operator, and F(V) — 0(|V| 2 ). We now prove that the spectrum of the operator 

9 1? (l + J B 2 )- 1 ^ f A 

near A = is approximately given by the linear dispersion curve Ai(£) with the associated eigenmodes given 
approximately by the Fourier modes exp(— iM/k). This follows from 

C l {e M,k ) = fc[£(e iM/fc ic^!) + (fcD(2flg«o + i(£/k)d^u ) - ic p d#u )e m / k ] , 

Of 

and therefore (p c mi J CAe lM/k ))(£) = x( — )iA[A',(0) - i£(2kDd&d k Vi + d#u ) + 0(£ 2 )]e m / k . Since 1 + B 2 (£) = 

I 1 

1 + <D(£) as a multiplier and d$e iM / k = i(£/k)e M / k we find 



J(d#(l + B 2 )- 1 p^ i e m/k )(£) = x(^)(AiW + 0(t))e M ' k . 



For notational convenience we diagonalize the linear part of (3. 30), (3. 31 ) by setting 



5" 



1 o\ hp' 
Si 1 U, 



(3.33) 



(3.34) 



where Si € C co ([-k/2,k/2),L(C,H m (T 27V ))) is a multiplier with suppSi C {h/8 < \£\ < 4/4} . Thus, v c = ip 
and P s v s — v s , and, by definition, 

5- 1 A5 = diag(A c ,A s ), (3.35) 
with X c (l) = x{f )(AiW + C(^ 3 )), cf. (|3.33|). In these coordinates, (|3.32| becomes 



d 4 « c = A c « c + cV™fA^> s ), 
a t w s = AV+P£ f Jv> c ,u s ), 

where A/" is a smooth nonlinear map from H™ +2 x H™ +2 into for every m > 1. 



(3.36a) 
(3.36b) 
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3.5 The moving frame 



To prove Theorems [T] and [2] we want to set up renormalization processes based on (3.36). For this we need to 
remove the 0(£) terms in Ai(£) = i(c p — c g )£ — a£ 2 + 0(£ 3 ). Therefore we define 



J- 1 



.u s 



^i(c p — c g )U 



(<d,e,t). 



This yields 



d t u c =\ g (£)u c + {d#+i£/k)p c ml rt(u c ,ii s ), 
d t u s =A s g (£)u s + P^AT(u c ,u s ), 



(3.37) 



(3.38a) 
(3.38b) 



where X g (£) = \i(£)— i(c p — c g )£ and A s g (£) — A(£) — i(c p — c g )£. The factors e ±lcH drop out of the nonlinearities 
since as multipliers they commute with the mode filters and 



(u* 2 )(£)e- icU = {{e icU v)* 2 ){£)e- iclt 
and similar for higher power convolutions. 



Kelt 



(v* 2 )(£), 



In general, (3.37) does not correspond to a simple transform in i?-space. However, if u has the special form 



u(t,l,ti) = a(£,t)g($) then, cf. (3.9), 



,fc/2 

v(d, t)= c ie ^/ k+ ( c "- c ^&(t, i)g(<&) d£ = a{tf/k + (c p - c g )t, t)g(ti) 

J-k/2 

= [a(x- c g t,t) +d#a(x- c g t,t)<j>^('&,t) + h.o.t]p(i9). (3.39) 
Thus, (3.37) will be responsible for recovering the group speed in Theorems [l] and [2j which motivates the index 



g in (3.38). On the other hand, completely transforming (3.36) to a comoving frame would make the linear part 



spatially and temporally periodic, and the subsequent analysis would require Floquet theory in time and thus 
be more complicated. 



The key features of (3.38) are the following. By construction, 

X g (£) = (A!(£)-i(c p -c g K) = -a£ 2 + 0(£ 3 ). 

We have 

{dv+'d/k^Nitf + u s ) = fj(£)N~ c (u c ,u s ), 



(3.40) 



(3.41) 



where \fj(£)\ = C£ and J\f c maps H m+2 (n) x H m+2 (n) into H s (n) for all s E N. In particular, by the calculations 
from |DSSS09j . 



f}(£)N c {u c ,u s ) = f3i£(u c 



h.o.t, 



(3.42) 



with p = — ^uj"(k), and where the higher order terms h.o.t are discussed later. The spectrum of A* is left of 
Kez < — (To, hence u s is linearly exponentially damped. Thus, heuristically, if for now we ignore u s and h.o.t in 
(3.42), then, as explained in S2.4 we have the following situations: in Theorem [l] and in Theorem [2] case (i) (with 
10" = 0), corresponding to Proposition 2.4 cases (i) and (ii), respectively, the whole nonlinearity is irrelevant and 



we obtain Gaussian diffusive behavior of u c ; for case (ii) of Theorem [2] (ui" ^ 0), corresponding to Proposition 
2.4 case (iii), the dynamics are governed by the Burgers equation for u c . 



The (unavoidable) drawbacks of the coordinates (3.38) are their relatively complicated derivation, and that (3.38) 



is quasi- linear while the original system (1.2) is semi-linear. 
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4 The results in Bloch wave space 



To prove Theorems [T] and [2j in |5]we set up renormalization processes for (3.381 in Bloch space. For this we 
need Bloch spaces with regularity and weights in £. Thus we first collect a number of definitions and basic 
properties. We recall that H m2 (nii)={u 6 L 2 (R) : ||n||ij-m 2 ( TOl )<oo} with ||it||fl™2( mi ) = ||up mi || ffm2 (R), where 
p(x) = (1 + l^l 2 ) 1 / 2 , and that F is an isomorphism between H m2 (mi) and H mi (m 2 ), where the notation 
H mi (m 2 ) = H mi (m 2 ) is used to indicate functions that live in Fourier space. 

Similarly, for L > and mj , m 2 , b > define 

B^{m 2 ,b) := {v E H^((-Lk/2,Lk/2),H^((0,2^))) : ||t>|| B ™i (ma , t) < oo}, 

II^IIb™ 1 ^^) = X! \\( d l d i} V )p b \\L 2 ((-Lk/2:Lk/2),L 2 (T2„))- 
Q<mi /3<m2 

Here p = p(£), i.e., we introduce a weight in the Bloch wave number £, and the subscript L indicates that the 
Bloch wave number varies in [—kL/2, kL/2]. For fixed L > the weight p is irrelevant since, due to the bounded 
wave number domain, all norms || • Wb™ 1 (m-y^x) an d II ■ lle™ 1 (m 2 ,b 2 ) are ec L u i va l en t, but the constants depend on 
b\, b 2 and L. The purpose of the weights is to take advantage of the "derivative structure" of the nonlinearity in 
the equation for u c , see (3.42), and Lemma 5.2 below. 

Let B mi (m 2 ,b) := B™ 1 (m 2 ,b). Then J is an isomorphism between H m2 (mi) and B mi (m 2 ,b), with arbitrary 
b > 0, see, e.g., jSch98bl Lemma 5.4]. We define the scaling operators 

Tl 1/L : £ mi (m a , b) -> S™ 1 (ma, 6), [n 1/L i]{$, I) = C(i?, i/L). 

Only £ is rescaled, and i? is not, and similar to ( 3.37[ ) this does in general not correspond to a simple rescaling 
of v. However, note that u c = u c (£,t) does not depend on ??, i.e., for u° Bloch space is identified with Fourier 
space, and in this case we have 

J-\n ljL u) = F-^TZ^lu) = LK L u, (4.1) 
i.e., concentration at £ = in Bloch space corresponds to spreading in Finally, 

\\Ki/ L v\\BZH2, b ) < CL b+1 ' 2 \\v\\ Bmi(2 ^ (4.2) 
and, for Bl ll (m 2 ,0) with m 1 ,m 2 > 1/2 and I € {-L/2, L/2), 



n 1/L {K L u*K L v)(£,x) = J 



1/2 

u{£ — Lm, x)v(Lm, x) dm 

! ■> 

L/2 

u(£ — m, x)v(m, x) Am =: L _1 (ii *£ v)(£. x). (4-3) 

-L/2 



This will be used to express the rescaled nonlinear terms, where henceforth we will drop the subscript £ in *l. 

To recall the heuristics, as a model for Theorems [T] and (in which the nonlinearities are completely irrelevant), 
consider the Fourier transformed version of dtu = ad 2 u, Ut—\ = uq, i.e., dtu = —a£ 2 u, which is solved by 
u{£, t) = Q-^( t ^ 1 ) al u{£, 1). Then, for any c G C, or more specifically c € K since we consider real valued functions 
u, f c (£) = ce~ al is a fixed point of the renormalization map 

et^fiHe-^ 1 - 1 ^)^ (4.4) 

Moreover, for L > 1 being sufficiently large, this line of fixed points is attractive in H 2 (2). To see this, write 
u(£) = f c {£)+g(£) with g(0) = 0. Then, using \g(£)\ < (£/L)\\dtg\\ c o (by the mean value theorem) and H 2 C 1 
we obtain 

I|e-^ 2(1 - 1/L2) ^i/X5l| 2 ff2( a) < CL-TgWn^). (4.5) 
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Thus, v(£,t) = u(t l l 2 l,t) — > / c (£) as t —> oo is the expected scaling for u c in Theorem [2|i) . Theorem [ijii) is 



also based on (4.4) but we have a nonlinear correction to the asymptotic profile as explained in £2.4 



Similarly, for any c G R, <? c (^) = ic£e Q ^ is a fixed point of the renormalization map 



g(» :i4e -< , (i-i/i , ) I5 I 1 



(4.6) 



and again this line of fixed points is attractive in H 3 (2) n Xo, where Xo consists of functions with zero mean. 
For this write u(£) = g c (£) + h(£) with d e h(0) = and use \h{l/L)\ < (£/L) 2 \\dfh\\ c o and H 3 ^ C 2 . Thus, 
v(£,t) = t 1 / 2 u(t~ 1 / 2 £ 1 t) — > <? c (£) as t — > oo is the expected scaling for u c in Theorem jl] The need for u G C 2 
also explains the higher weight in x in Theorem [T] 

Theorem 3 [Diffusive stability]. Let Uo(-',k) be a spectrally stable wave-train and b G (0,1/2). There exist 
e,C > such that if \\{u c ,u s )\t=i\\H 3 (2)xB 3 (2,2) ^ £ i an d u c (0, 1) = ^Fk J u c ($, l)di? = 7 then the solution 
(u c ,u s ) to (3.38) exists for all t > 1, and there exists a ipum G K such that 



\\t^u*(t-V 2 £,t) - iTh im £e-"* 2 \\ H3{2) < CrW, 



(4.7) 
(4.8) 



Theorem 4 [Diffusive mixing of phases]. Let uo(-;k) be a spectrally stable wave-train and b € (0,1/2). 
There exist e, C > suc/i £/iai /or |^| < e i/ie following holds. 

(i) Assume that (3 = —^u>"(k) = 0, \\u c (£, 1)\\h 2 (2) < £ with w c (0, 1) = (f>d/{2nk), \\u s (-, l)||e 2 (2,2) < £ 
P s u ;5 (-,1) = u s (-,l). Then the solution (u c ,u s ) to (3.38) exists for all t > 1, and 



||^ 17 ^ c (£, t) - u c M\\h^2) < cr l / 2 +\ 
\\n lM u s (eMB*<2,2)<ct-v 2 + b , 



(4.9) 
(4.10) 



where u% (£) = 4>de at . 

(ii) If (3 = —^u>"(k) 7^ then the same result holds with u%{£) replaced by 



u%{£) = T 



ia ze 



13 1 + zed [&/Vka) 



(4.11) 



where ln(l + z) = ^<fid 



Before proving Theorems [3] and [4] we show that they imply Theorems [T] and [2] 
Proof of Theorem [TJ Given initial data in the form ( 2.15[ ) from Theorem [TJ i.e., 

u(x,t)\ t=0 = u (9 - 6» + (j>o(z); k) + v (x) with ||<fo||ff3(3)) ||uo||j?3(2) < £> 

we first need to extract (u c ,u s )\t=i and show that they fulfill the assumptions of Theorem[3j Then we translate 
(4.71,(4.8) back into (<f>,v) coordinates. 



Thus, as explained in Remark 3.2 let u(x, i)|t=i = tto(#; k{\ + d#<f>a)) + iuo($) with 9 = i3 — 4>o(fi) and 

woi.fi) ■= uoifi; k) - u (-d; k(l + d^fo)) + v (x). 



W.l.o.g. assume that (1 — P s )wq = 0, otherwise redefine <fto = P^f^o- This fixes the non-uniqueness in (2.15). 
Also, (j>o G H m (3) for all m G N due to the compact support of (f>o, and, with ipo = d^(po, J from (3.12) and S 

iu c ,u s )\t=i = JS- 1 i^ ,w ) = J^ Q ,w ~ SiVo) (4.12) 



from (3.34), 
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is well defined and fulfills ||(fi c ,w s )|| H 3 (2)xe 3 (2i2) < C x e and u c (0, 1) = ^ / u c (i?, 1) M = 0. 

1 i? 



We now use (4. 7), (4. 8) to recover Theorem 



Using (3.10) and -FfVe )W = - 



/47ra 2a 



- fl2 /( 4 «) we have 



where 



lim 



t1l t i/2U c {d, t) - ^tie-® 2 '^ 2 ^ = T- 1 t 1 / 2r R t -i/2U c {t, t) - f0 lim £e- Q£2 {&), 



: , and from Ci||ti|| ^(j^) < \\u\\ Hm ^ n) < c 2 \\u\\ H ™( m ) we obtain 



Vina 2afc' 

||ift tl/2M c (tf,i) - ^m0e-* 2 /< 4fc2Q > 11*2(3) < Ct-^ 2+b . 

Then, with 

i>{'d,t) ^u c (d + k(c p - c g )t,t) and tu(#,t) = 5iu c (i?/fc + (c p - c 5 )t, i) + w s (t?, t) 
we obtain, in L°° , 



r" 5 2P?/h- rv /■*+fe(c P -c g )t 



-i 2 /(4k 2 at) 



y/t J-oo V 2ak2f 

i.e., 4> rim = -4fcVa 3 7nAiim- Also tt;(i?,t) = since suppSi C {£i/8 < \£\ < £i/4} and 



(4.13) 



|K(<M)|| 



j M ' k u s {dJ,,t)dt 



-k/2 



r 1 ' 2 / ^"^/^(^r^tJfi+^a + O" 2 ^ 

-kt 1 / 2 /2 



(4.14) 



Finally, 



i? =6* - O + <t>{-&> t) = 6-B - 2fc 2 t- 1 / 2 V'ii m ae- ( ' 5+fc(c ^ c ' )t)2/(4Qfe2t) + 0(t~ l ) 
=0-9 o - 2k 2 t-^ 2 ^ lim ac- {x - c ^ 2 / < - 4at ^ + 0{t- r ) 

using = kx — wt = k(x — c p t) and the implicit function theorem. 



□ 



Proof of Theorem [2j First, assume that /3 = 0. As above we write u(x, t)|t=o = wo($; fc(l + <9^0o) + u>o($) 
with w (i?) := ito(i?; /c)— w,o($; fc(l + (9^0o)) + wo( a; )j an d where now ||0 o (')llff 2 (2) < £ an d 4>oW — > 0± as i9 — > ±oo. 
Again, w.l.o.g. assume that (1 — P s )wq — 0. Then, for 

we obtain 2irku c (0, 1) = J u c (i), 1) di9 = <fid, and Theorem [4] applies to 

({t c ,u s )| 4= i = J<S _1 (V>o,wo) = J(^a,w a - SiV>o)< 
Thus, with ^>($, £) = u c ($ + fc(c p — c g )t, and £) = 6>iu c ($ + k(c p — c g )t, t) + u s ($, t) we obtain, in 



\/4irk 2 t 



l5+fc(Cp— Cg)t 



+ (0+ - 0_)erf((x - c g t)/yfAat) + O^ 1 ) 

and «;(#,*) = ©(t^ 1 / 2 ^) as in | |4.13| and j4.14| ) above. Hence 

= + = + </-_ + (0+ - (p-)erf(x-Cgt)/V4ai) + 0{r 1/2 ) 

and by shifting the 0(i _1 / 2 )-part to v we obtain part (i) in Theorem [2] Part (ii) with (3 ^ works in the same 
way. □ 
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5 Renormalization 



We first prove Theorem [3j the minor modifications needed to prove Theorem [4](i) are then explained in ^ 5.5 
while the changes for the slightly more complicated proof of Theorem Qii) are explained in [5.6 



5.1 The rescaled systems 

Based on the asserted behavior t 1 / 2 u c (t~ 1 / 2 l,t) — > itpn m £e~ ae we introduce, for n € N and L > 1 chosen 
sufficiently large below, the variables 



Then {u c n ,u s n ) fulfill 



where 



(5.1) 
(5.2) 



(5.3a) 
(5.3b) 



u°(k,t) := L n u c ( K /L n ,L 2n T) = L n [n L - n u c ](n, L 2n r), 
u s n {ti,K, r) := L n u 3 ($, n/L n , L 2n r) — L n [R L - n ii s ](-&, k, L 2n T). 

d T u c n {K,T) - ~X gin (K)u c n (K,T) = L 3 ™A/^«,0(k,t), 
d T u s n (ti, k, t) - A g , n u s n (ti, k, t) = L 3n Af°{u c n ,u s n )tf, k, r), 

A ff ,„(/s) = L 2n \ g ( K /L n ), A 9) „ = L 2n K L -„k g K L ^ 

K«,<) = fj(K/L n )n L - n M%L- n n Ln u c n ,L- n n L nu s n ), 
K(u c n ,u s n ) = n L - n P^M(L- n n Ln u^L- n n L ,u s n )- 

Except for the different scaling due to m°(0,^) = 0, (|5.3[) has a very similar structure as, e.g., |Sch96[ eq.(30)] or 



|Uec041 eq.(3.2)]. Thus, similar to (2.381, we shall consider the following iteration: 



solve (5.3) for r e I := [L~ , 1] with initial data " \ (■&, k,L~ 2 ) =L[ r ' 

\ui \ui 



(5.4) 



Formally, (5.3) is solved by the variation of constant formula, i.e. 



+ / e ( T -^-WL 3 x c «,<)( K)S )d S , 

Jl/L 2 



- p (r-l/L 2 )A s K t » f , S 



K,T) 



1/L 2 



9 ^L Nn{u c n , u s n ){&, K, s) di 



(5.5a) 



(5.5b) 



However, (5.5) can not be used to construct the solution since (5.3) is a quasi-linear system, as it can be seen 
from J\f : H m2 (mi) x H rn2 (mi) — > H m2 ~ 2 {m,\) in (3.36). To solve (5.3) we use maximal regularity methods 



[LM6S] for parabolic equations in (weighted) Sobolev spaces as in Uec07|. A posteriori, (5.5| can then be used 
to estimate the solutions. Thus, we first note some properties of the linear semigroups and the nonlinearities in 



(5.5), and then explain how to obtain local existence for (5.3) 



5.2 Estimates on the linear semigroups and the nonlinearities 



We shall need some detailed estimates on the linear semigroups and the nonlinear terms in (5.5). The idea is 



to exploit the derivative-like structure in the Bloch wave number k of Af c as expressed in (3.42 ) by relaxing the 
weight, and to regain the weight using e^ r_T ) A s-™. Thus, from this point on, the weights in k become important. 
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Lemma 5.1 There exists a C > such that for all L > 1 we have 

|| e (r-rOA g ,„~c ||s ^ (2 ^ < Cmax{lj ( T _ TO" 672 } K || B 3„ (2,2-6), (5.6) 

l|e (r - r ' )A -<|| si „( 2 , 2) <Cmax{l, (r - r')- m2/2 }e-^ L2 "( T - T ')||<|U3 (2 _ m2 , 2) . (5.7) 



Proof. Equation (5.6 1 holds since the real part of A 9 .„(k) = L A s (k/L™) = — an +0{k ) is bounded from above 



by the parabola — aoK 2 , while (5.7) holds since A gi „ is a relatively bounded perturbation of L 2n (d$ + m/L n ) 



and by construction has spectrum left of — L 2rir y . □ 

The following lemma transfers the fact that derivatives give higher powers of upon rescaling to general 
convolution operators with a "derivative-like" structure. 

Lemma 5.2 Let mi E N, 7 > ; and K G C™ 1 ([-1/2, 1/2) 2 , ff 2 (T 27r )) iwi/i ||^(/c-^, £)|| ff 2 (T2ir) < C(|k-^| + |€|) 7 . 
TTien 

,L/ 2 

(v» H> (M l/L K) (v,w)(k, ■&) := / [Tii/jr Jf]( K - £4 -!?){;(«;, - d^ 

defines a bilinear mapping (Mi/lK) : B™ 1 (2, 2) x (2, 2) /3™ 1 (2, 2), and i/iere exists a C > suc/i iftai for 
all L > 1 we raaue 

||(Xl/L^)(«»|| er(2 , 2 _ 7) < C7i- min tT ' 1 > 1 1 « I Ib" 1 (2,2) 1 1 ^ 1 1 B2* 1 (2,2) " 



Proof. This holds due to sup £ 



PL ' 



□ 



(i+e 2 )i/ 2 

Lemma 5.3 Let \\{un, Un)\\[B 3 Ln (2,2)] 2 ^ Rn < 1- There exists a C > smc/i t/iat 

L 3 »||A£«X)|| B 3 (24) < CL-"i? 2 . (5.8) 
TTie term L 3n Af„ can be split according to the number of $ derivatives as L 3n Af^ = + A/"* 1 + Af^ 2 such that 

\\K,i\\B3(2-i,2) < CR 2 n . (5.9) 



Proof. We write L 3n JV,1 = si + s 2 , where, as explained in { 3.5 the lowest order term si in L 3 ™A/^(m° , u s n ) reads 



Sl («) = L 3n i/3 — ft L -„(iT"7^<r 2 ( K ), 



(5.10) 



cf. (3.42). This yields || «i || b 3 (2.1) < CL ™i?„ by direct calculation. The remaining terms s 2 can be estimated in 



a similar way using Lemma 5.2 and taking into account the finite support of (d^+i£ / typ^N \u c + u s ) in Fourier 
space. 



This does not work for L 3n J\f^(u^, u s n ). However, here we do not need an additional factor L n , and (5.9 1 simply 
follows by checking the number of derivatives in Af and using (4.3 ). □ 



5.3 Local existence 



Since (3.36) and hence (5.3) is quasilinear we cannot combine Lemma 5.1 and Lemma 5.3 to directly show local 



existence for (5.3 1 via (5.5 ). Instead we use maximal regularity theory from |LM68) . For / = (r , T\) and r, s > 
let 

H r > s {I, mi) = L 2 {I, H r ( mi )) n H S {I, L 2 (mi). 



Since (3.36 1 is a parabolic problem these spaces only occur with s = r/2 and we set K m2 (I, mi) = H" l2 ' m2 ^ 2 (mi). 
Then, for any given weight b > 0, Bloch transform is an isomorphism between K m2 (I,mi) and 



K mi (I, 772 2 , b) = L 2 (I, B mi (m 2 , b)) n H m * /2 (I, B mi (0, b)) 
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Similarly, for every n, let 

AT™ 1 (I, ma, b) := n 1/L ,K m ^ (I, m 2 ,b) := L 2 (7, (m 2 , 6)) n i? m2 / 2 (I, B™ 1 (0, 6)), 



i.e., the subscript L™ again indicates that the Bloch wave number varies in [—kL n /2,kL n /2). From (4.2) we 
have 



(/,m 2 ,6) 



(5.11) 



Recall that for each n the weight 6 in k gives an equivalent norm in K™^ (7, m2, 6), but the constants depend on 
n. We also need subspaces of functions that vanish sufficiently fast at To, and define 

K m2 {I, mi ) := {v G K m2 (I, mi ) : d J T v(-,T ) = for j eN,j < m 2 /2 - 1/2}, 
J?p(/,m 2 ,6) := G K™i(I,m 2 ,b) : -,r ) = for j G N, j < m 2 /2 - 1/2} 



We set 



I=(iT 2 ,l), 



and for (m^, u^)| T=L -2 G (2, 2)] 2 construct solutions (m£,m*) G [Ar£„ (/, 3, 2)] 2 to ( 5.3 1. Note again that for 
u c n we can identify Bloch space with Fourier space such that in fact u c n G AT 3 (7,2) (in the Fourier sense) with 



supp-S^(r) c I n = {\k\ < L n £i/4}. We abbreviate (5.3) as C n U n =M n (U n ), where 

c jj = ( d T V,l(K,r) -A g , n («)<(«, t) \ 

\ d T u s n r) - A fl) „u*(i?,K,r)J ' 



(5.12) 



and, for m 2 > 2, we first consider the linear inhomogeneous version of (5.3) with zero initial data, i.e 



C n U n (T) = JV„(t), Af n G [ if£»(/,m 2 - 2,2)] 2 , C/ n | T=i - 2 = 0, (5.13) 
where moreover for the first component Af£ of N n = (J\^,J\^) we assume 

A> G AT 3 (7, 2) (in the Fourier sense), and suppA>(r) C I n = {\k\ < L n £x/A}. (5.14) 

Lemma 5.4 There exists a C > 0, independent of n G N , smc/i that for all Af n G [o AT£„ (J, m 2 — 2,2)] 2 which 
fulfill \5. 14 ) there exists a unique solution of (5.13) with 



WUnWa^ (i,m 2 ,2) - C\\Jvn\\ Kl n (I,m 2 -2,2)- 



(5.15) 



Proof. The first component 9 t m^(k,t) — A 9! „(k)m^(k, t) = A/"^ is independent of $ and thus can be solved by 



the variation of constant formula using (5.6 1 (with 6 = 0). For the second component we use resolvent estimates 
for the solution of 

(A-A 9 , n )<=A^. 



There exists a C > such that for m 2 > 2, A^ G B™£ (m 2 - 2, b) all A G C with ReA > we have 



i rillB™ 1 (m 2 ,6) 



\M w ' /2 \\u\\ Kim < c(|KI| srn x (m2 _ 2jb) + |A| (m - 2) ||A>||^ ;(0 , b) ) (5. 



16) 



Similar to Lemma 



5.1 



this holds since A ffi „ is a relatively bounded perturbation of L (d$ + in/L n ) and by 



construction has spectrum left of —L 2n ~f . See also |Uec07[ Appendix A. 2] for an explanation of how to obtain 



resolvent estimates in weighted spaces. From (5.16) we obtain (5.15) by continuation of J\f£ for r G K, Laplace 
transform, and the Paley-Wiener Theorem. In fact, in (5.16) we could choose A to the right of — L 2 ™7o, but 
ReA > is enough to show (5.15) with C independent of n. □ 
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We denote the solution operator of (5.131 by qC^. To solve the nonlinear problem we write U n = V n + W n 



where V n € 3, 2) is a continuation of U n \ T=L -2, which exists due to [LM681 Thm 4.2.3]. Then W n fulfills 

C n W n = G n (W n ), W n \ T=L -2=0, where G n (W n ) = N n {V n + W n ) - C n V n . (5.17) 

5.4 and 



The idea is to show that for W n € qK^„(I,3, 2) we have G(W n ) € qK^„{I,1,2) and use Lemma 
estimates on the nonlinearity to apply the contraction mapping theorem to 



We set 



<s>(w) := c- l G n (w n ). 



Pn ■= \\(&n>Un)\T=l\\[Bl n (2,2)]= 



(5.18) 



(5.19) 



and obtain the following local existence result, taking into account that (tt£, u s n )(1/L 2 ) and (w^_ 1 ,uf l _ 1 ) are 



related by {u c n ,u s n )\ T=L -2 = LKi/ L (u c n _ 1 ,u s n _ 1 )\ T= \ and hence, by ( [42] ), 

||«,<)| T=i - 2 || B 3 (2,2) <^ 7/2 Pn-l- 



Lemma 5.5 There exist Ci,C 2 > 0, independent of n, such that the following holds. If p n —\ < C\L 7 ^ 2 ; then 

(5.20) 



there exists a unique solution (u^jUfJ € [i^£„ (I, 3, 2)] 2 to (5.3) with 

\\{K,K)\\{Kl n (i,3,2)p < C 2 L 7/2 p n - 1 . 
Moreover, for all t\ > L~ 2 and any wi2 G N there exists a C3, independent of n, such that 

ll«X)ll[KS(( ri ,l), mai 3)]» ^ C3i 7/2 Pn-l- 



(5.21) 



Proof. From standard Sobolev embeddings we have that A/" n is a smooth mapping from K\ n (I 7 3,2) to 



K^ n (1, 1) 2), see also Lemma 5.3 To show that G n (W n ) in (5.17) is in oi^£n(/-l,2) we have to fulfill one 



compatibility condition, namely G n (W n )\ T =L- 2 = 0, which holds by construction. For sufficiently small p n -i, $ 



is a contraction since M n is quadratic and higher order. In particular, combining Lemma 5.4 with a slight adaption 



of (5.9) to the time dependent case we find that C\, C2 may be chosen independent of n. The higher regularity 
follows by a standard bootstrapping argument: for almost all r s (L 2 , 1) we have (u£,u*)(r) G S£„(3, 2). 



Starting again at such a r the required compatibility conditions to apply Lemma 5.4 are automatically fulfilled. 
This yields fl5.21|). □ 



5.4 Proof of Theorem [3] (Diffusive stability) 



Due to the loss of L 7 / 2 in Theorem 5.5 we need to improve |5.20 ) to iterate (5.4). Given a local solution , u 



with the higher regularity (5.21 ) this will be achieved by using the variation of constant formula and a suitable 
splitting of u c n . 

For u c € H 3 (2) we define Uu c — d K u c (0), which, by Sobolev embedding, gives a continuous map, i.e. |IKt c | < 



C||it c || fj3, 2 \- Here, and also in (5.25 ) below, we need the smoothness in the Bloch wave number, which for u° we 



again identify with the Fourier wave number. To prove Theorem [3] we write 



';(«, 1) = iip n g(K) + r c n {n), u s n ( K , d, 1) = <(??, k), 



(5.22) 
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where g(n) = ne aK and r£(0) = d K r^(0) = 0. This makes sense since u n (0,r) = for all n E N and all 
t G [1/L 2 , 1] if u c (0, 1) = 0. Substituting (|5.22| into (JKB yields 



< = e {1 - L ~^ n Ln 1/L r n - X +I c n + Res„, 

< = ^- L -^ K 'LHy L r n -l + I S n , + I S n,l + I S n,2, 

where, using the notation L 3n Af£ = J\f£ + M£ i + A/"^ 2 from Lemma 

e (i-)A 9 ,^c ( -c )<)(r)dT) T . j= r 



(5.23a) 
(5.23b) 
(5.23c) 



5.3 



H = l 



1/L 2 



1/L 2 



s (X-r)A,. B ^. ( <)<)(r)dT) 



and where the residual in (5.23d) is defined by Res n = i^ n _ie^ L ^ Xg " LlZx/ L g — h\) n g. We also define 

Pn,c = IKIU?„ (2,2) and Pn,s = lknlUf„ (2,2)i 



which gives, cf. (5.19), p„ = ||{£ || B s (2)2 ) + IKIU^ (2,2) < C|V>. 



n| 1 Hn,c 1 Pn,s 



Now assume that p n _i < L 7 / 2 . Then from (5.61,(5.8) we immediately obtain 

\2 



\lpn ~ V>«-l| < CL- n {C 2 L 7 l 2 p n ^Y 



(5.24) 



with C 2 from (5.20). Moreover, \\e^~ L )x '- n LK 1/L g - g\\ B 3 n{2 . 2) < CL~ 2n and hence ||Res n || S 3 nC2 ,2) < 
CL- 2 " |t/;„_i I . Next, we have 



l|e (1 - L " a)X "-i^ 1/i r=_ 1 || fl 3 (2 , 2) < CL-'Wr^l 



,(2,2)- 



(5.25) 



This follows from r < f l _ 1 (n/L) = (^) 2 9 2 t^_ 1 (k) for some k between and k. Here again we need the smoothness 
in k. Combining the above estimates we arrive at 

Pn.c < CL- X p n _^ c + CL- n (C 2 L 7 ' 2 Pn ^f + CL- 2n \i> n ^\. (5.26) 



To estimate p n>s first note that 



,(2,2) 



<Ce- yoL2 " J L 7 / 2 p n _ 1:S <L-V«-i,. 



for L sufficiently large. Next, 1^ and i"® l can be estimated using (5.7) and (5.9) to 



ll J n,ollsi„(2,2) + H-fn,llls3„(2,2) < CL "(^Vi-l)" 



(5.27) 



(5.28) 



-I In 

However, for the quasi-linear part 2 we have to use the higher regularity (5.21 1 and split J£ )2 = fij L 2 •" dr + 
Jl/2 ' ' ' dr to obtain 

r l/2 r l 

II^, 2 |Ib£„ (2 , 2 ) <C(C 2 L 7/2 Pn-i) 2 / (l-rl-V^'^'dr+C^Vi) 2 / e^^dr 

L Jl/L 2 Jl/2 

<C{C 2 + Ct)L- n {L 7 l 2 p n ^f. (5.29) 

This yields 

Pn,s < CL- X p n _ x + CL- n (L<l*p n _ x f (5.30) 
Now, let L > Lq with Lq so large that CL^ 1 < L~^ b ' and let po — \\(u c ,u s )\\b^( 2 , 2 ) < £ -4 . Then, combining 



(5.24), (5.26), (5.27) and (5.30), iteration shows that there exists a ipiim S K such that 

\lp\lm ~ 1pn\ + Pn,c + Pn,s < L~" (1 ~ b) as 71 -> OO, 



(5.31) 



where the correction L nb takes care of the powers C n arising in the iteration. This discrete convergence implies 
Theorem[| using t = L 2n T and the local existence Theorem |5.5[ □ 
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5.5 Proof of Theorem |4](i) (Diffusive mixing, Gaussian case) 

The main difference compared to the proof of Theorem [3] are different scalings for u°,'U*, which are now based 



on (4.4 1 instead of (4.6). Thus, we introduce 



i£(K, t) = Tl L - n u c {K, L 2n r), u' (0, K, t) = K L -„u s (i9, k, L 2n r). 



We want to show that 



We obtain 



u° n {K, 1) -> d e 



I* (k, 1) — > as n — > oo. 



d T U c n [K,T) - \ g:n {K)u c n {n,T) = L 2n Af°[u c n ,u s n )(K,T), 

d T U^(-d,K,T)-A g<n U S n (-d,K,T) = L 2n Af r l(u c n ,{i s n )(l},K,T), 

where A Si „(/«) = L 2n \ g (n/ L 2n ) and A S) „ = L 2n TZ L - n A g TZLn as before, and 

K(u c n ,<) = f)( K /L n )K L - n tf(n L nul,K Ln u s n ), 



(5.32) 



(5.33) 



(5.34a) 
(5.34b) 



The renormalization process reads 



solve (5.34) for r £ I := [L 2 , 1] with initial data 



t=L- 2 



5.5 



(5.35) 

=i 

with «,<)(i/L 2 ) e 



The local existence for (5.34) works exactly as for (5.3). In contrast to Lemma 
[B£„(2,2)] 2 it suffices here to take (u° n , u s n )(l/I?) £ [5f„(2,2)] 2 in order to extract the asymptotics (|5.33|) 



cf. (5.45). Thus, we set 



Pn ■- IK"^'", S l )|r=l||[s2„(2,2)]2, 

and for p n _ 1 < CiL~ h / 2 wc obtain a local solution (u c n ,u s n ) £ [K 2 n (J, 3, 2)] 2 to (5.34) with 

\\«,0\\[iCl n (I,3,2)]2 < C 2 L 5/2 Pn - 1 , 
and for each t\ > L~ 2 and m 2 £ N there exists a C3 such that we the higher regularity 

ll(<X)ll[i^ n (( n ,l),m 2 ,2)P ^ C 3 £ 5/2 /Jn-i- 
The estimates for the linear semigroups work as before, i.e., here 



| e (r-r')A s ,» fi c|| (22) < Cmax{lj ( T _ T 



-^-M- b /2}|| 



u nllfl?„(2,2-b); 



| e (r-r )A„n fi -|| (2i2) <Cmax{l, (r - r 



The nonlinearities are now estimated as follows. 



/\-m 2 /2i -7 L ™(t-t') ||„~.s II 



(5.36) 

(5.37) 

(5.38) 

(5.39) 
(5.40) 



Lemma 5.6 Lei || (w^, ) || [bJ^ (2,2)] 2 ^ Rn < 1. There exists a C > smc/i i/iai 

L 2 "|KK,<)|U= <CL-»i? 2 . 



(5.41) 



The term L 2n J\f^ can be split according to the number of d derivatives as L 2n J\f^ = M£ + 1 + ■M? 2 such that 



\\KiWBl^-i,2)<CR 2 n . 



(5.42) 
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Proof. Apart from the different power counting, the proof works like the one of Lemma 5.3 with the crucial 
difference that now the term s\ from (5.101 vanishes since j3 = by assumption. □ 



Similar to (5.22 1 we now set 



u c n (K,l) = <t> d g{K) + r c n (K), u s u {k,<&,1) = r£(«,0)i 



(5.43) 



where g(n) = e aK and r£(0) = 0. Here no variables 4> n are necessary since, due to the conservation of total 
phase shift, i.e., d t u c {0,t) = for all t. We obtain 



'Ki/lCi +I c n + Res n , 



r s — n^-L 2 )^i'T? , T r s -A- T s -A- T s 4- T s 



(5.44a) 
(5.44b) 



where Res„ = </>d(e^ L 2 ^ Xg - n 7li/ L g — g) and 



It = L 



l/L 2 



e (l-)A 9 ,^c ( -C !<)(T)dT) j. j= / e (MA 3 ,^ jKi<)(T)dTi 



l/L 2 



with L 2n Afn =J^n,o + J ^n,i + K\.2 from Lemma 5.6 Clearly, ||Res„|| B 2 n(2]2) < C\(/) d \L- 2n , and 



IS 2 (2,2) 



(5.45) 



which follows by writing r^_ 1 («/L) = (^)9 K r^_ 1 («;) for some k between and k. Combining this with (5.391 

Pn,c < CL- x p n _ ltC + CL- n (C 2 L 5 / 2 Pn ^) 2 + C\cf> d \L- 2n , 



and (5.41) thus yields 



and similarly 



Pn,s < £ _1 P»-i,a + C(Cl + Cf)L-"(L 5 / 2 fti -i) 2 



(5.46) 



(5.47) 



The proof of Theorem |4jji) now follows by iteration. At this point the assumption \(j)d\ < e is crucial. This can 
be seen by computing p n:C in powers of L for n = 0, 1, 2, . . . starting with po.c = 0. Hence, we need \4>d\ < L~ d 
with d sufficiently large. □ 



5.6 Proof of Theorem |4](ii) (Diffusive Mixing, Burgers' case) 

Essentially, Theorem Qii) is again based on the scaling (5.32). However, the crucial difference to the case j3 = 
is that now the analog of ( |5.41 ) no longer holds. Therefore, we need to scale u s differently, i.e., we blow up in 
order to avoid problems with the quadratic terms involving u s in the critical part, in which the term i£[3(u c * u c ) 
will give the Burgers dynamics for u c . Thus, for small p > we introduce 



££(«,t) =K L - n u c ( K ,L 2n T) 1 k,t) = L n(1 -P)K L - n u s (ti, k, L 2,1 t), 



(5.48) 



to obtain 



d T U r n (K,T) ~ X g:n (K)u C n (K,T) = L 2n Af^(u C n ,U S n ){K,T), 

d T u s n (ti, k, t) - A g , n u s n (#, «, r) = L n( - s -rttf*(u c n , k, r), 
where again A 9i „(/c) = L 2n \ g {nj ' L n ) and A g ,„ = L 2n lZ L -n AgTZ^ , but now 

K(u c n ,u s n ) = ^(«/L")^ i -jO'(7i i n< J i-»( 1 -rt^ i »«*) ) 



(5.49a) 
(5.49b) 
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Accordingly, the renormalization process reads 



solve (5.49) on r G I := [L , 1] with initial data 



. ft? 



=L- 2 



(5.50) 



The estimates for the linear semigroups are again (5.39) and (5.40), while the nonlinear terms are estimated as 
follows. 

Lemma 5.7 Let \\(u^,Un)\\[B 2 „(2,2)] 2 < Tt n < 1- There exists aC > such that L 2ra A/^(w° , ) = S1+S2+S3+S4 
wz£/i si = if3n(v,n * u£) and 



II s 2||b£„(2,1) < CL ll«nllB2„(2,2)' 
IMU| n (2,l) ^ C ' i_n(1_P) ||W™llsi„(2,2)IKIIsi„(2,2), 
IMU M 2,1) < CL- 2 ^ 1 -")^. 



(5.51) 



The term L^'P^Af* can be split according to the number of D derivatives as L™*- 3 p 'A/" T f = A/* + A/^ j + A/"* 2 
such that 



T»,illB?„(2-i,2) 



'(1-P)|| 7 --,C||2_ 1 



R n)- 



(5.52) 



Proof. The term S2 contains the quadratic terms in u c n except for * i.e., si is of the form S2(k, t) = 

h(n/ L n )(un * u^) with /i(k) = 0(k 2 ). S3 contains the quadratic interaction of u c n and u s n , and S4 contains the 



remaining terms. Then (5.51 ) follows from the finite support of A/j£ in Fourier space. It is in S3, S4 that the blowup 



scaling uffi, k,t) = L n ^- p ^U L - n u s ('&, k, L 2n r) is useful. Equation (5.52) again follows by straightforward power 
counting. □ 



The terms involving only in (5.52) blow up as n — > 00. However, combining (5.52) with the exponential 



damping in the stable part, cf. (5.40), we still get local existence for (5.49) with constants independent of n. We 
let 



Pn '■= ll(Wn)"n)|r=l||[B= (2,2)P) 



(5.53) 



and for p n -i < C\L 5 / 2 obtain a local solution {u^ L ,u s n ) € [^, 2 (I, 3, 2)] 2 to (5.49) with ||(Un)^n)ll[K' a n (7 3 2)] 2 — 
C 2 L 5/2 p n -i, as in (5.371, which moreover enjoys the higher regularity \\(ul_,u s n )\\ [k 2 n{{Tul) m2 2 - )]2 < C 3 L 5 / 2 p n _ 1 , 



cf. (5.38) 



Similar to (5.22) and (5.43) we now separate from u c n the lowest order asymptotics, now obtained from the 



Burgers equation. However, due to the contribution s\ = \(3k,{u'^ 1 * u^) of the nonlinearity to the asymptotics 
here we work out an intermediate step and split u c n in a r dependent way. In detail, let 



where 



"n( K > T ) = «n,*( K > T ) + «„(/S,r) 
«n,»(«> r ) =X(K>/L n )uZ(K,T) with U%(t,t)=J : 



(5.54) 



/? 1 + zerf (V/y/kai) 



0. 



cf. (4.11 ), and x from (3.17). Consequently a n (0, r)=0 for all n, r due to the conservation of total phase. Then 



d T a c n (K,r) - \ g , n (K)a c n (K,T) = L 2n {Af^(u c n ,u s n )( K ,T) - A^K^, 0)) +Res„, 
d T uW,K,T)-A a , n ZM0,K,T) = L"( 3 -P)A>«,<)(^, K ,r), 



(5.55a) 
(5.55b) 



where 



Res r , 
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Lemma 5.8 There exists a C > such that sup £ -2< r<1 ||Res n (r)|| B 2 n ( 2 ,2) — "I'M- 
Proof. By construction, i.e., since u c n t is an exact solution of the Burgers equation, 

Res„Kr) = CL- n (o(K 3 )u^ +0(k 2 (u c ^ * u c n J)^j 

which can be estimated in B\ n (2, 2) by L" n \4>d\ since u%(n, r) is analytic and exponentially decaying. □ 
Now setting 

K( K > 1 ) = K,*( K ' 1 ) +r n( K ): 1) =<(«,!?), (5.56) 

the remainder of the proof of Theorem |4|ii) works as the proof of Theorem |4ji) in j ]5.5| □ 
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